Asymmetrically Doped Polyacetylene 
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Doped one- dimensional (ID) conjugated polymers, such as polyacetylene, have a conductivity of 
some metals, like copper. We show that when this polymer is asymmetrically doped, unexpected 
properties are revealed, when compared to the behavior of the symmetrically standard doped systems 
(SDS). Depending on the level of imbalance between the chemical potentials of the two involved 
fermionic species, the polymer can be converted into a ID partially or fully spin polarized organic 
conductor. 
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I. INTRODUCTION 



The interest in polyacetylene was revived recently by 
many reasons. From the theoretical point of view, it 
can been considered a ID version of graphene, at least 
with respect to the fractionalization of the electric charge 
PiSQ- From the experimental and technological sides, 
the advent of the .emzconrfuctor nanotechnology i with 
new possibilities of developing semiconductor devices 
made from polyacetylene Q strongly motivated the study 
of transport properties of conducting ID polymers [6|. 
Polyacetylene is a ID chain of CH groups with alternat- 
ing single and double bonds. The system is half filled, 
i.e., has one electron per site and, in the tight binding 
approximation, it would be a metal. However, due to the 
interaction of the electrons with the lattice or, the sponta- 
neous Peicrls dimcrization, the polyacetylene is an insula- 
tor. After doping, the conductivity of trans-polyacetylenc 
is greatly increased, reaching metal- like properties [7|. 

The model which describes the electron-phonon in- 
teractions is the celebrated Su-Schrieffer-Heeger (SSH) 
Hamiltonian whose continuum version is known as 
the Takayama-Lin-Liu-Maki (TLM) model [l]. The 
TLM model is a relativistic field theory with two-flavor 
Dirac fermions. We shall use the Gross-Ncveu (GN) 
model [lol | in 1-1-1 dimensions since its Lagrangian is 
equivalent to that of the TLM model [ll|, and will pro- 
vide us with analytic solutions. In the field theory con- 
text, the discrete chiral symmetry is spontaneously bro- 
ken in the GN model (at zero temperature and chemi- 
cal potential) , and there is the dynamical generation of a 
mass (gap) for the electrons. The (global) external chem- 
ical potential ji is introduced in the theory to represent 
the extra electrons that arc supplied in the system by the 
doping process. At a critical chemical potential = 

where Aq (« 0.7eV for trans-(CH)x) is the constant 
band-gap, the GN model undergoes a first order phase 
transition to a symmetry restored (zero gap) phase p^ . 
The critical doping concentration [lj| yc can be related 
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with the critical chemical potential as yc = api^ [l^ . 
where iV (= 2) is the number of spin degrees of free- 
dom of the (delocalized) tt electron, Ti is the Planck's 
constant divided by 27r, vp = kpfi/ra is the Fermi veloc- 
ity, kp is the Fermi wavcnumbcr, m = /2toa?', and a 
(= 1.22 A) is the lattice (equilibrium) spacing between the 
X coordinates of successive CH radicals in the undimer- 
ized structure. io(~ 3 eV) is the intercarbon transfer 
integral for tt electrons, or simply the "hopping param- 
eter", resulting in vp ~ lO^m/s for trans-(CH)x, which 
is of the order of the velocity of the Dirac fermions in 
graphene [lB|- Employing the GN model in the large- 
N (mean-field) approximation, previous work [3| found 
a very good agreement with the experimentally found 
y^ (= 6%). Recent calculations including finite correc- 
tions to N improved the theoretical prediction for y^ 
and showed that j/c is solely weakly affected by thermal 
effects [11. 



The novelty is brought about by an "asymmetrical 
doping" , defined as an imbalance between the chemical 
potentials of the electrons with the two possible spin ori- 
entations ("up" =t, and "down" =[) inserted in the sys- 
tem by the doping process. Since the densities of the t 
and J, electrons are directly proportional to their chem- 
ical potentials, we investigate the consequences of the 
asymmetrical doping in the magnetic properties of the 
frans-polyacetylene. 

Experimentally, the chemical potential asymmetry can 
be achieved by the actuation of an external magnetic 
field on the system, which breaks the spin-1/2 SU(2) 
symmetry. We show that with low imbalanced doping 
(with respect to the critical doping), the energy gap is 
zero, although there is an imbalance between the ] and 
I electron populations. The situation changes substan- 
tially with strong imbalanced doping, favoring one of the 
movers from the two spin orientations (f or |). In this 
case, the system acquires a small but stable non-zero gap 
Ao((5/^) < Aq, and there is an absence of the not favored 
spin oriented electron. We argue that this may have ob- 
servable consequences in the magnetic properties of the 
asymmetrically doped polyacetylene (ADP). 
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II. MODEL LAGRANGIAN 

The Lagrangian density of the TLM model in the adi- 
abatic approximation (i.e., neglecting the lattice vibra- 
tions) is given by 



N 



^TLM = ^i'^ (idt - ihvp^sdx - loM^))^' (1) 
' -A^(x) , 



2nhvFXTLM 



where tp is a two component Dirac spinor tJj^ 



representing the "left-moving" and "right-moving" elec- 
trons close to their Fermi energy, respectively, and j is 
an internal symmetry index (spin) that determines the 
effective degeneracy of the fermions. We define 1 =t, 
and 2 =J,. The gamma matrices are given in terms of 
the Pauli matrices, as 70 = cri, and 75 = — (T3. A(x) is a 
(real) gap related to lattice vibrations, Atlm = ^^j^ is a 
dimensionless coupling, where a is the 7r-elcctron-phonon 
coupling constant of the original SSH Hamiltonian, and 
K is the elastic chain deformation constant. 

The equivalence between the TLM and the Gross- 
Neveu (GN) model, is established by setting Atlm = 
^S^. The partition function of the GN model in the 
imaginary time formalism [17] is given by 



Z ^ [D4'''][Dijj]exp / dr dx {L 



^GNJ 



(2) 



where P = 1//cbT, ks is the Boltzmann constant, and 
Lgn is the Euclidean GN Lagrangian. The (CH)x can 
be doped after synthesis both by using chemical or elec- 
trochemical techniques. If an electrochemical process is 
being used to add electrons through external circuits, 
the polymer is being reduced. When the doping re- 
moves electrons it is said that the polymer is being ox- 
idized. In the standard way, the chemical potential is 
introduced in the theory by adding to I/gn the term 
fiip^ipi which means that the same quantity of t and | 
electrons is being inserted in the system. The asym- 
metrical doping is established by adding to Xgn a term 
/ij-i/iT^'i/jT _|_ jj^ijji'' ipi- ^ with ji^ ~ p, + Sfi, and = p, — Sji. 
We will also choose p, ^ fic since, as we will see next, 
this is the chemical potential at which the minimum of 
14//(A, /i|,l) jumps (when Sfi = 0) from Aq to p^ . 
This allows us to write Lgn = + + Lq, where 

= ^ta {-lodr + ihvF-fA - A + joni,i) V't.i, ^0 = 

~ 2nvp\GN ^'^ ■■ ^ = "^^^0 = ■0^cri, and 71 = ia2. Integrat- 
ing over the fermion fields we obtain 



"f^gn" ; m^^detD 



'1 ' 



where Di,2 = -i/3[(-iw„ + ^|_|) - -uf7o7iP - 7oA], and 
LOn = (2n -I- l)7rr are the Matsubara frequencies for 



fermions. Summing over n and carrying out the deter- 
minantal operation, we find 



InZ = -13V 



1 



-A' 



V 



+00 



dp 

2TTh 



+I3E+ + 2\n[l + e^'^'^T j + (3E:f + 2 In [l + e^'"'^ 
-I3E+ + 2 In (l + e^^^^t\ + + 2 In (l + e^^^^ 



where V is the ID "volume", E'^ ^ = Ep±^-\,i, and Ep = 
\Jv\p^ + A^. The "effective potential" for a constant A 
field is given by 14// = ^■^^ InZ, then 



(3) 



-A' 



V 

+°° dp 

2Trh 
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fesTln (1 + e"'^^^) + keT In (l + e~'^^T" 
keTln (1 + erf^^t) + fesTln (l + e^'^^r 



(4) 



The zero temperature and chemical potentials Ve// is 
written as 



VeffiA) 



1 



UvfXgn 



A^- I % 2E 



(5) 



where we have introduced a momentum cutoff A to reg- 
ulate the vacuum (divergent) part of T4//(A). It is easy 
to see that Sfi has lifted the degeneracy of the conduction 
{+) and valence (— ) bands in the matter part of Vgff. 
The renormalized effective potential turns out to be 



^.n(A|, <0) 



where mp is an arbitrary renormalization scale. The 
minimization of 14// (A) with respect to A gives the well- 
known result 



Aq = mpe^ 



(7) 



At finite chemical potentials and in the zero tempera- 
ture limit, we have 



K//(A,/iTa) 



hvpX 



' -A^-r^2E„ 



GN 



(8) 



T^fi Jo 



where ph^ = 1 — A^ is the Fermi momentum of 



the moving electron. We integrate in p, observing 
that the renormalization is the same as before, to obtain 
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Veff (A,MT,i) 



27rhvF 



AHn 



^GN 
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- A2 



mtV^t 



27rhvF 



A^ln 



We see that the minimum of Vef fhyp/ Aq for < Sfic 
is at A = 0, while for Sfi > 6fic, a narrow gap A(J/x) 
(9) "appears" through a first order (quantum) phase tran- 
sition. However, this transition now is from A = to 
A((5/i) ^ 0, i.e., is in opposite direction from that ob- 
served in SDS ((5/i = 0), as a function of /i. The behavior 
of Vef fhvp/ Aq shows that the gap will eventually vanish 
for very high asymmetry, or (5/x >> S^c- 



- A2 



Mi 



In 



Ao 




sion: 



Ao(Ai) = 6(/ic - fJ-)AQ. 



(11) 



A^ - 2^l^AlA 



At^O. 



(12) 



where 8i_2 = ©(m^ x~ A^) is the step function, defined as 
Q{x) ~ 0, for a; < 0, and Q{x) = 1, for x > 0. Minimizing 
Ve//(A, fJ-i^i) with respect to A yields the trivial solution 
(A = 0), and 



(10) 



where we have made use of Eq. ([7]) to eliminate m^? in the 
gap equation. The ground state is determined by jointly 
finding the solutions of the equation above Ao((5/i), and 
the analysis of the effective potential at the minimum, 
Fe//(Ao(5/x)). In the SDS, defined for Sfi = 0, the solu- 
tion A = Aq represents a minimum of Ve//(A,/Lt) while 
fi < ^c, (where the critical chemical potential, /ic, is ob- 
tained by Veff{A = 0,/ic) = Veff{A = Aq, /ic)), and 
A = is a local maximum. When fi > fic, A{fi > fic) is 
a local maximum and A = is turned into a minimum 
through a first order phase transition [l^, [3] , agreeing 
with experiment [l^. Thus, the gap as a function of the 
chemical potential in the SDS has the following expres- 



The situation changes significantly in an asymmetri- 
cally doped system, where A = still represents the min- 
imum of Veff{A, /Xf^i) for (5/i < (S/ic = 0.38Ao, where S^c 
is obtained from the equality 14// (A = 0, ^^(J^c)) = 
Ve//(A = Ao((5^c),MT,i(<5/^c)), with Ao((5^c) being the 
solution of Eq. (fTO|) for ^'^^^{Sfic)- For S^/ Aq > Sfic, the 
02 function prevents the "J," term in the effective poten- 
tial of participating in the minimum. Thus, Eq. (jlOp can 
be rewritten as 



We show in Fig. ([T]) the (non-dimensional) function 
Veffhvp/Al for (5/i/Ao = 0, Sh/Aq = 0.3, S^/Aa = 0.5, 
(5/^/Ao — 0.6, and 5^/ Aq = 0.7, as a function of A/Aq. 




A/A„ 



FIG. 1: The top curve is the effective potential of Eq. (9) 
for <5^/Ao = 0, and jj, — fic- The second curve from top to 
bottom is for Sjj,/ Aq — 0.3, with the minimum of 14// at A = 
0. The following curves are for Sfi/ Ao — 0.5, Sfi/ Aq — 0.6 
and (5/i/A() = 0.7. There are stable gaps now at Ao(5^i = 
0.5Ao) > Ao{SfJ. = O.6A0) > Ao(5^ = 0.7Ao) all of them 
agreeing with the solution of Eq. (|12|l . 



III. MAGNETIC PROPERTIES 

As we mentioned, the imbalance in the chemical poten- 
tials of the t and | electrons can be caused by the applica- 
tion of a static magnetic field B = Bok, perpendicular to 
the wire, with a Zeeman energy AE = SzgnsBo p^.l20t. 
where Sz = ±1/2, g « 2 is the effective g— factor, 
fiB = eh/2m « 5.788 x 10"^ cV T'^ is the Bohr magne- 
ton, m is the bare electron mass, and Bq is the magnetic 
field strength. In the present case Sfi — \AE\ = fisBo. 

The number densities n^^^ = — ^ 14// (A, /i-f,j,) are 
obviously imbalanced due to the asymmetry between 
and Hi, and will depend on Sfi. Before the critical asym- 
metric doping, Ao((5/i) = 0, and the densities read 



1 



Trhvp 



(13) 



For such a low imbalance, compared to the critical chem- 
ical potential asymmetry Sfj^c, the total number density, 
rix, is the same as in the symmetric limit or, in other 
words, is independent of the applied field: 
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fJ-c 



(14) 



This partial spin polarization results in a net magnetiza- 
tion 1201 of the chain 



that clearly docs not have the simple form of the Pauli 
susceptibility. 

Finally wc find the critical magnetic field separating 
the partially and fully polarized magnetic phases Sfic = 
0.38Ao = mbSo,c, yielding 



M {Sn < dfj,c) 



nhvp 



s^l (15) 



Bn = Mn. 



The magnetic susceptibility in the low asymmetrical dop- 
ing regime is 



0,c 



4.6 kT, 



(22) 



which is a magnetic field of very high intensity, compared 
to the maximum current laboratory values, that are 
up to 20 T (DC-niagnetic fields), and up to 331 T 
(mega-gauss fields) j21| . 



dM _ 2fj.l 



OBq irfivp^ 



(16) 



which is exactly the Pauli expression for noninteracting 
electrons. 

Increasing the asymmetry beyond the critical value 
means that, effectively, only |-moving electrons of the 
conduction band are present in the system. In terms of 
this "higher" chemical potential asymmetry the densities 



n|((5^ > (5^c) 



ttHvf 



(17) 



where A ~ A{Bq) in the equation above is the solution 
of Eq. for a given /i| , and 



niidfi > Sfic) = 0, 



(18) 



meaning that the system is fully polarized. The total 
number density is now 



nriSfi > Sfic) = n| = — — V(/ic + fisBoY - IS?. 
nnvF 

(19) 

In the fully spin polarization regime the magnetization 
of the chain reads 



IV. SUMMARY 

To summarize, we have investigated the mean- field 
zero-temperature phase diagram of a ID polymer upon 
asymmetric doping. Our main finding is that an ADP 
can be converted into a partially or fully magnetized or- 
ganic wire, depending on the level of doping imbalance, 
i.e., the difference between the t and | electron densi- 
ties, induced by an external magnetic field. Although the 
(critical) magnetic field i?o,c necessary for a fully magne- 
tization of the polyacetylene be of very high magnitude, 
we have seen that any non-zero Bq produces a net mag- 
netization. This doable partial spin polarization can be 
a manifestation of the interesting "itinerant magnetism" , 
for which conduction electrons are responsible. It would 
also be desirable to study the polarization mechanism 
in the trans- (CII)x at finite temperature, and the subse- 
quent calculation of thermal effects on the magnetization, 
critical magnetic field, and susceptibility, within the field 
theory approach. We hope to address these issues in fu- 
ture. 
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M{Sfi > Sfic) = fJ-Bn-\ 



sj{fic + fiBB^y- - A2 > Afo. 



The magnetic susceptibility in this regime is given by a 
rather involved expression 



xiSfi > Sfic) 



^^B 



(Mc + f-LBBa)fii 



OA 
'dB„ 



TTflVF ^(^t, + flBBo)^ - A2 
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